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Abstract 

A family of systems related to a linear and bilinear evolution of roots of polynomials 
in the complex plane is introduced. Restricted to the line, the evolution induces 
dynamics of the Coulomb charges (or point vortices) in external potentials, while its 
fixed points correspond to equilibriums of charges in the plane. The construction 
reveals a direct connection with the theories of the Calogero-Moser systems and Lie- 
algebraic differential operators. A study of the equilibrium configurations amounts in 
a construction (bilinear hypergeometric equation) for which the classical orthogonal 
and the Adler-Moser polynomials represent some particular cases. 
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I. Introduction 

In the present paper we propose to discuss a Bilinear Hypergeometric Operator 

H£[f, g] = (f"g - 2A/V + AV7) P + \ (f'g + A 2 g'f) P' + (f'g - Ag'f) U + Xfg 

P ■= p( z ) = A + Bz + Cz 2 , U :=U(z) =a + bz, AeK (1) 

f:=f(z),g:=g(z), f := ^ 
and to study integrable dynamics: 

of roots Xj, yj of polynomials in a complex variable z 

m n 
Qm(z, t) = Y[(Z - Xi(t)), p n (z, t) = JJ(^ - Vi(t)) 
i=l i=l 

induced by the action of (1) 
alq m dp r , 



(2) 



Pr, 



Aq m -^ = [p n , q m ], X nm = (Am - n) (V + (n - km)P" /2) (3) 



alt alt 



The above construction has a nice physical interpretation: The fixed points of (2) correspond 
to equilibrium distributions of n and m Coulomb charges (or point vortices in hydrodynamics) 
of values 1 and —A respectively in external potentials on the plane or cylinder, while the real 
solutions describe their motion on the line or circle. 

It should, however, be noted that the physical analogies are not complete, because alx*/alt and 
aly*/alt should appear in place of alxi/alt and alyi/alt in the lhs of (2). Nevertheless, all equilibrium 
solutions in the plane as well as time dependent solutions on the real line are same for (2) and 
corresponding physical systems 

The electrostatic interpretation for roots goes back to works by Stieltjes on the classical or- 
thogonal [20], and by Bartman on the Adler-Moser polynomials [4]. These results represent special 
cases in our construction. 

The Bilinear Hypergeometric Equation 

H£[f,g](z) = (4) 
is a natural extension of the Gauss hypergeometric equation 

H£[f, l](z) = P(z)f"(z) + (U(z) + \p'(z)^ f + Xf(z) = 



and the recurrent relation 

9 i+1 ](z) = 9': +1 (z)9 t (z) - 29^)9'^) + 9 i+1 {z)9'l{z) = 0, P = 1,U = 

for the Adler-Moser polynomials 9i(z) [1]. Another special case 

A = 1, p( z ) = - z 2 , U = 

of (4) provides an interpretation for Huygensian polynomials in two variables studied by Y.Berest 
and the author in connection with the Hadamard problem [6]. 
The paper is organized as follows: 

In the next two sections we study a linear evolution of polynomials, which is a particular case 
of the bilinear dynamics. We show that such an evolution corresponds to dynamics of pairwise 
interacting charges iff it is induced by second order Lie-algebraic (hypergeometric or "quasi exactly 
solvable" ) differential operators. The system of charges then can be embedded into an integrable 
Hamiltonian (in general, elliptic Calogero-Moser or Inozemtsev) model related to a Coxeter root 
system. Skipping consideration of the "quasi" and elliptic cases in the sequel , we classify remaining 
cases by types of corresponding hypergeometric equations. 

Sections IV- VI are devoted to the study of the bilinear evolution and its fixed points in the spe- 
cial case A = 1: Introducing the bilinear hypergeometric operator, we show that its action induces 
dynamics, which, in some settings, can be interpreted as a motion of unit positive and negative 
Coulomb charges (or a system of vortices in hydrodynamics) in external potentials. We find that 
it can be embedded in a flow generated by a sum of two independent Calogero(Sutherland)-Moser 
Hamiltonians. In this way (2) can be integrated by the Lax method if A = 1. 

Considering polynomial solutions to H\ [f, g] — in section V, we analyze equilibrium con- 
figurations of charges. It turns out that such solutions can be obtained from associated linear 
problems by a finite number of the Darboux transformations. 

In section VI we discuss degenerate limits of the bilinear dynamics related with the Kadomtsev- 
Petviashvilli equation and give an interpretation to a set of algebraic solutions of a particular type 
of (4) obtained earlier in connection with the Hadamard problem in two dimensions. 

In section VII we introduce polylinear (/-linear) hypergeometric operators and related dynami- 
cal systems, with (l)-(3) presenting a special case 1 = 2. In this picture, / distinct types of charges 
move in external potentials, interacting with each other. Such a dynamics can be embedded in 
a Hamiltonian system of I species of particles of a Calogero-Moser type. In contrast with the 
I — 2, A = 1 case, there is no separation of the Hamiltonian flow in independent components and 
the Calogero(Sutherland)-Moser type potentials are not related to the Coxeter root systems. 

In section VIII we return to generic two component dynamics, presenting some arguments in 
favor of integrability of (2) for arbitrary real A. 

Some open questions are discussed in the concluding section of the paper. 

II. Linear Evolution 

This section, which is a generalization of works by Choodnovsky & Choodnovsky [10] and Calogero 
[8], is devoted to the study of a particular case of the bilinear evolution. 
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Let 

V = Span{M,.--,* n_1 ,3 n } (5) 

be a linear space of polynomials over the complex numbers C, V = C n+1 of the degree less or 
equal to n in z. 

We consider the evolution of polynomials p(z, t) 

= L\p(z,t)}, p(z,t)=T(t)f[(z-z i (t))eV (6) 

i=i 

under the action of a time independent linear operator L e End(V, V). 

Rewriting (6) in terms of roots Zi(t),i = 1 . . . n and the common factor T(t) we arrive to the 
following 



Lemma 1 The linear evolution equation (6) is equivalent to the following dynamical system 
T- 1 ±T = T(z 1 ...z n ), (a) 

(7) 



z i = v{z i \z 1 ,z 2 ,...,z i ,...,z n _ 1 ,z n ), i = l...n, (b) 



v is a rational function symmetric in the last n — 1 (all but Zi, the hat in (7a) denotes omission 
of Zi) variables. 

Proof: Representing L in the matrix form 

n 

L\z l \ = Y,U,z\ L tj eC (8) 

3=0 

we equate the lhs and the rhs coefficients at different powers of z in (6). Equation (7a) is obtained 
by picking out a coefficient at z° = 1. Expressing dT/dt from (7a) and substituting it to the rest 
of equations we get: 

= fi(zi...z n ), i = l...n (9) 

where fi are polynomials symmetric in z iy i = l..n and <7j stand for elementary symmetric polyno- 
mials 00(24 ...z n ) = l, a 1 {z 1 ...z n ) = X;™ =1 Zi, 0-2(21 ...z n ) = J2 i<:j ZiZ j: . . . 

Equation (9) is a linear system for dzi/dt with determinant nKj^ ~ z i) • ^ ^ s no ^ singular, 
provided all distinct, and (9) has a unique solution. 

Since fi are symmetric in all arguments, this completes the proof. 

We call (7) a system with two body interactions if 

v(z i \z 1 ,Z 2 , ■ ■ ■ ,Zi, . . . ,Z n -i,Z n ) = ^ w(Zi,Zj) +u(zi) (10) 

In (10) and in the sequel the following notations iV := {1, 2, ...n}, M := {1, 2, ...to} are used. 
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Theorem 1 : A system generated by (6) for n > 2 is a system with the two-body interaction if 
and only if L is a (modulo adding a constant) second order differential operator 



L = P(z)£ + U(z)£ - *U'(z) - ^P"(z) 

P{z) = A + Bz + Cz 2 + Dz 3 + Ez 4 , U(z) = a + bz + cz 2 - 2(n - l)Ez 3 

with polynomial coefficients P(z) and U (z) at most degree four and three respectively. 
Under condition (11), (7) becomes 



(11) 



A 7 1 

i=-^E^-»w d2) 



In other words, (12) is the most general system with two body interactions induced by a linear 
evolution in the polynomial space. 

Remark: The cases n = 1, 2 are excluded from the theorem, since any linear operator induces 
a two-body dynamics. 

Proof: Introducing the following set of differential operators 

Cij = J^f[( Z Tz- k )^ £>*\A = ^z\ Q<i )3) k<n 

we see from (8), that any L G End(V,V) can be represented by a differential operator with 
polynomial coefficients Qi(z) 

n n n ,j 

i = EE L ^ = E^)^ ( 13 ) 

i=0 j=0 i=0 

Substituting (13) to (6), using (7a) and imposing condition (10) we have: 
Q := t(z!, z 2 , . . .) - Qo(z) - Y.i (Qi( z ) ~ u i z i)) U 

-2 Ei<j (Q2(z) - w(z h Zj)) Utj - 6 Ei< j<fc Q^z)Utjt k - = 

where tj = l/(z — zA. 

Taking the n + 1 th partial derivative dzdzid d Z2 q z — = we get the following equation 



{t 1 ...t n ) 2 UQ n {z)Y,U-Q' n {z) \ =0 

V ieN J 



Picking out the coefficients at different powers of t^i G N, we find that Q n (z) = 0, since z,ti....t n 
is a set of independent variables. 
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Proceeding by induction and taking the n, n — 1, ...4th derivative of Q 

d n Q 



dzdz ix dz i2 . . . dzi 



0, i\ < %2 < ... < ij, 2 < j < n 



we eliminate all Qi with % > 2. 

Hence , the operator L is at most of the second order. 
Decomposing L into the sum of homogeneous components 

L = Y L u L { : Cz j -> Cz i+j , L { = a iZ l+2 ^ + hz i+1 ^- + 

ct2;^ dz 

i 

and remembering that V is spanned by z l ,i = 0..n only, we have 

L l [z n - 2 ) = L t [z n - 1 ] = L l [z n }=0, for i>2 

It is immediate that 

a>i — bi = Ci = 0, for % > 2 
Proceeding in this way, we impose similar conditions for % — 1 and i = 2 

I^z"" 1 ] = L 2 [z n ] = L 1 [z n ] = 

getting the most general expression (11) for a second order differential operator L e End(V, V). 

The sufficient condition of the theorem is proved by the direct calculation: Substituting (12) 
and 

d n 

T = L n n + L n _i n Zi + L n _ 2 ,n z i z j 

i=l i<j 

into (11), expressing the matrix elements in terms of a,b,c,A,B,C,D,E,n and using the 
following identity 

i 1/1 i \ 

(14) 



we show that (6) holds identically. This completes the proof. 

Although (12) is not a Hamiltonian system, the following proposition holds 
Proposition 1 (12) is a trajectory of a system with the Hamiltonian 

^ = 2£^.^g^ + (^ 



(15) 



In other words, the Hamiltonian equations of motion 

'd 2 Zi\ P'(zi) (dzi\ 2 dV 



P{ Zi ) 2 \ dt 2 J 2P{zi) 2 \dt J dz, ' 

are corollaries of (12). 



i e N (16) 
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Remark: Equations (16) are reduced to the Newtonian form 

d 2 ^ dV 



dt 2 dfa' 

by the change of variables 



% E N (17) 



* = «*>• ^r=vk (18) 



Proof of proposition 1: Expressing second derivatives through (12) = 2~2jeN ~lk~§r (IfO ^ 
direct calculations, we evaluate the lhs of (16): 



dzi \ . ^ (zj ^ Zi- z 3 2 ^ P(zi) 



Using the identity 



£ -*)(*„-«) = 5 + ( " " 2> \ E{Z ' + Z ' ,+ D(Z ' + Z ' ) + 3 C J 



we get (15), which competes the proof. 



III. System Classification, Lie- Algebraic Operators and 
Calogero— Moser Models 

Consider the n + 1 dimensional representation of the Lie algebra sl(2, C) 
by differential operators 

J+ 

acting on (5). 

Operator (11) is an element of the universal enveloping algebra of s/(2,C), L = ^ a vt; J v + 
j3< ; J <: ,'v,S = ±,0 (19). Such Lie-algebraic operators are called "quasi-exactly solvable" [21]. They 
can be separated in nine nontrivial equivalence classes under the linear-fractional transformations 
of the independent variable z. Based on the invariant-theoretic classification of canonical forms 
for quartic polynomials [12], operators (11) can be placed in the nine canonical forms with 

(1) P(z) = l (4,5) P( z ) = 1tz 2 

(2) P(z) = z (6,7) P{z) = {l T z 2 ) 2 (20) 

(3) P(z) = -z 2 (8,9) P(z) = (z 4 + tz 2 =f 1) 



[J°,J±] =±J ± , [J-,J+] =2J° 



2 d n d n d 

z—-nz, J — z- -, J = — 

dz dz 2 dz 



(19) 
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P(z) 


z = z((f>) 
(18) 


Root 
System 


Calogero-Moser 

-I — v , . -ITT/ / 1 \ \ -\ /-. ■— \ 

Potential V(z(<f>)), (15), (17) 


Polynomial Eigen- 
functions of L (11) 


1 


Z = (f) 




E,< fc (^7F + |E,(« + ^') 


if b 7^ Hermite 

#n (* + f )) 

if 6 = 0, 
l,z 


z 


z — 2_ 

Z 4 


BC/D„ 


4,4 

^j<fc (<t>k-<t>i) 2 Ofc+<M 2 

, V 2a 2 , ^ 
^ ^ 8 


if 6 7^ Laguerre 

n V / 

if6 = 0,a= 1-j 


-z 2 


z = exp(i<p) 




V • , ,j 1 + 1 V • ae~ i<t>j + 6e*^' 


z n L b-2n+l ^ 


1-z 2 


z = cos(0) 


BC/D„ 


V 1 I 1 

Z^j<fc sm 2 (4> k -<t>j)/2 sm\4> k +<t>j)/2 
, (6+a) 2 ab 
' j 2 sin 2 4>j 2 cos 2 ^ /2 


Jacobi 

/ a+b+2 a-2-b\ 
Pn 2 ' 2 \z) 



Figure 1: Four generic classes of Hypergeometric systems 



Most general Hamiltonian systems (15) in this classification are elliptic Inozemtsev models (for 
trigonometric and rational Inozemtsev Modes see e.g. [15], [19] ) related to A n and BC/D n Coxeter 
root systems. 

Skipping analysis of "quasi" and elliptic cases in the sequel, we content ourselves with "exactly 
solvable" hypergeometric operators dealing with the first five classes only (and linear U (z) = a+bz) 
in (20). 

Figure 1 gathers needed information on the first four (the fifth one is a hyperbolic version of 
the fourth) classes, which are related to the rational/trigonometric Calogero(Sutherland)-Moser 
systems. 



IV. Bilinear Evolution, A = 1 

In this section we introduce a special case H\[-,-] := of the Bilinear Hypergeometric 

Operator (1) and study related integrable dynamics of roots. 
Let 

Vi = c n+1 , v 2 = c m+1 , v 3 = c m+n 

be linear spaces of polynomials of degree less or equal to n, m and n + m — 1 respectively. 
Consider the evolution 

q^-p^ = H x [ P ,q}, P eV 1 ,qeV 2 (21) 
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under the action of a bilinear operator H\ : V\ x V 2 — > V3 on the monic polynomials 

p=H(z-x t (t)), q = l[(z - Vi (t)) (22) 

ieN ieM 

of the nth and mth degrees. 

Lemma 2 The bilinear evolution equation (21) is equivalent to the following dynamical system 

j^Xi = vi(x i \x 1 ..x i ..x n \y 1 ..y m ), i = l...n 

f t Vi = v 2 (yi\yi..y i ..y m \x 1 ..x n ), i = l...m 

where v\ and v 2 are rational functions (v\ is symmetric and y . v 2 is symmetric in 

y x ..yi..y m and x). 

Proof: is similar to the linear case (lemma 1), except that the polynomials are monic now. 
dx/dt and dy/dt are uniquely expressed from a linear system of equations with the determinant 

Ylj<ieN( x j- x i) YljeN,ieM( x j-Vi) Yli^eMiVi-Vj)- lt is 11011 singular provided all roots are distinct. 
Again, we study systems with two body interactions which means that 

vi(xi\x\y) = J2j^i w n( x i, x j) + J2ijWi 2 (xi,yj) 

v 2 (yi\y\x) = J2j^iW22(yi,yj) + 'E iJ w2i(y i ,x j ) 

As in the linear case, it is natural to look for differential operators inducing integrable dynamics 
in a system with two-body interactions. It turns out that such operators exist and are extensions 
of the linear case. 

Proposition 2 The bilinear operator H\ nm : V\ x V 2 — > V3 

Hx nm [ P , q] = (p"q - 2p'q' + pq")P + ±(jpfq + pq')P' + (p'q - q'p)U + \ nm pq (23) 

P(z) = A + Bz + Cz 2 , U(z)=a + bz, 
X nm =(m-n) (V + (n- m)P"/2) 

induces dynamics 

= ma) (- eu,^ ^7 + sr=i ^-) - - m*) 

(24) 

iVi = 2p (^) (E7= m* ^7 - E?=i ^7) - U( Vi ) + \P\yd 
by action (21) on (22). 
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Remark: Equation (21), (23) may be written in the form of the Schrodinger evolution equation 
with a time-dependent potential 

d A = P (zW + (U+ 1/2P'W + (P'(ln q)' - 2P(ln q)" + A) ^ (25) 

where ip — p/q. In this setting, we study the time evolution of polynomial q and a rational function 
tp with denominator q, which is a rather inconvenient formulation for our purposes. 

Remark: Substituting P(z) = 1, U(z) = —k, (dk/dz = 0), and A = in (25) and reexpressing 
it in the formally self-adjoint form we obtain the nonstationary Scrodinger equation 

^ = $"+W$,W = -2(logr)", (26) 
r = q 

which is a second equation of an auxiliary linear problem for Kadomtsev-Petviasvilly hierarchy 
[18] (with q as a r-function). The solution to (26) is now a quasirational function 

$ = E. exp(kz + k 2 t) 

We observe similarities with Krichever construction [16] , [17] for the rational Baker-Achieser 
function. In more details, the B acker- Akhieser function ^(z,t,k) is a special n = deg(p(z)) = 
m = deg(q(z)) case of the above quasirational function 

*=*=(i+i:^)-p(^) 

with divisor of simple poles defined at points x i: i E N, N = M. 

Proof of the proposition 2: Is essentially similar to the proof of sufficient condition of theorem 
1: one substitutes (24) in (23) and uses identity (14). 

Similarly to the linear case, equations of motion (24) can be expressed in the Newtonian 
coordinates (18) 

dd>i dH d6j dH , s _ , „ N . „ T . 

i<jeN ieN i£N,j£M i£M i<jeM 
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Should we have d(f>*/dt, d9*/dt instead of dfa/dt and dOi/dt in the lhs of the equations of 
motion, the system (27) would be a Hamiltonian system of n positive and m negative vortices or 
Coulomb charges on the plane or cylinder [2], [3]. It is not Hamiltonian in our case, but has the 
same fixed points in the plane or cylinder and dynamics on the real line or circle. 

Let us discuss the question of integrability of (24). 

In the linear case, which is a particular case m = of (21), there were two lines of approach 
to the integration of system (12): 

For the first possibility, the linear dynamics (6) in the finite basis (5) allowed us to find p(z, t) 
(and Zi(t)) at any t. 

The second way to find was to solve the Hamiltonian system (17) with initial conditions 
given by (12) itself. 

Obviously, the first of above approaches does not apply in the bilinear case, since the evolution is 
not linear any more. Therefore, we use the second method, trying to embed (24) into a Hamiltonian 
system. 

Lemma 3 If the odd function $(x) satisfies the functional equation 

$(x)$(y) + $(z)$(x) + $(y)$(z) = (29) 
whenever x + y + z = 0. Then the following identities hold 

T ^ x ) = H Yl '$2*(x n -x i )&(x n -x j )=2 H^-Xj) 2 (30) 

ieN,i=/=njeN,j=/=nneN i<jeN 

I 2 (x,y) = 2j2H ~ VmMxi - x j)- 2 J2Y1 ^ " x mMVi - Vj) 

meM ieN jeNj^i rneN ieM jeMj^i 

+ ^ _ X mM%i - Vj) - ®( X i ~ VmMVi - Xj) = (31) 

rneN ieN jeM meM ieM jeN 

Proof: is a calculation. 

Theorem 2 (24) can be embedded into the flow generated by the sum of two independent Hamil- 
tonians 

H = H + + H_, 

2 i / j \ 2 



(zi-zj) 2 2^ P{ Zi ) ' ±v ' v ' 2 



In other words , the Hamiltonian equations of motion 



d* dv + (cm d* _ dv.jge)) ( . 



where ( is given by (18), are corollaries of (24)- 



11 



Remark: Some results related to the special case P(z) — 1 of theorem 2 (rational A n root sys- 
tem) were obtained by Veselov [22] , who studied rational solutions of the Kadomtsev-Petviashvili 
equation. In particular, it was found that poles of (unbounded at infinity) rational solutions of 
the KP equation (which are coordinates Xi,i G N in our case) move under the Calogero-Moser 
flow with nonzero external potential. It is interesting to note that the nondegenerate external 
potentials U + and LL coincide only in the above mentioned special case. 



Proof of theorem 2: Let us check that the Hamiltonian equation of motion for x 

1 / d 2 X,:\ P'(Xj) ( dx,:\ 2 8 



P(xi) 2 \ dt 2 J 2P( Xi ) 2 V dt J dxi 
holds, expressing the second derivatives through (2) 



V + (34) 



d 2 Xi dxj d ( dx,i \ dyj d ( dxi 



dt 2 ^ dt dxj V dt J ^ dt dyj \ dt 

J= l 3=1 y] V 

By direct calculations we get 

lhs(34) = - A (2W 1 (A, B, C\x, y) + 4W 2 (A, B,C\x)+ li " < '' 



d Xl ^ — - — v-,-,-.-, 2 p{xi) 

where 

+ Y Y Y P<yVj ^ Y YY P<yXj ^ 1 YY 

w 2 (a,b,c M = e E E (z 

and = A + Bz + Cz 2 . 

Let us evaluate W^A, B, C\x, y) = AW^l, 0, 0\x, y) + BW^O, 1, 0\x, y) + CWi(0, 0, l\x, y). 
In W\{\, 0, 0|x, y) we immediately recognize identity (31) with <3>(x) = Consequently 

W^i(l,0,0|a;,y) = / 2 (a;, J /) = 

Changing variables = exp(0j), y^ = exp(#j) we find that 

Wi(0, 0, l|x, y) = J 2 (0, 0) + Anm(l - n + m) = Anm(l —n + m), = coth(x). 

Finally, using linearity of W\ with respect to parameters A,B,C we write W\(0,l,0\x,y) 
W^O, Q,i\ x +iy+±)- W^O, 0, l\x, y) - \W X (0, 0, l|x, y) = 
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Therefore 

W!(A,B,C\x,y) = Anm(\ - n + m) 

Applying (30), we evaluate W 2 in a similar way, getting (34) with V + given in (32). 
The proof is completed by applying a similar procedure to yi. 

Corollary 1 (1) is integrated by the Lax method 

Proof: Since (32) are Calogero(Sutherland)-Moser Hamiltonians, equations of motions (33) 
can be represented in the Lax form [19]. Thus solutions to (24) may be found from (32) subject 
to initial conditions given by (24) itself. 



V. Equilibrium Configurations, Bilinear Hypergeometric Equa- 
tion 

The fixed points dxi/dt = and dyi/dt = of (24) describe equilibrium of the unit positive and 
negative Coulomb charges in two dimensional electrostatic or point vortices in hydrodinamics re- 
spectively [20], [3]. The polynomials p and q (22) must then satisfy an ordinary bilinear differential 
equation 

H Xnm [p,q} = (35) 

which is a special case A = 1 of the Bilinear Hypergeometric Equation (4). 

Studying the dynamics of roots we supposed that they are distinct and polynomials p and q 
do not have common factors (lemmas 1, 2). However, solutions of (35) may have multiple roots 
or/and common factors. In this circumstances we need to modify (24). 

Remark: one does not encounter such a problem in the linear case since polynomial solutions 
of ordinary Hypergeometric equation (classical orthogonal polynomials) do not have multiple roots. 

Proposition 3 Let p and q be polynomials of orders n and m satisfying (35) and 

n m 

p/q = p/q, p = Y[(z-x i y i , q = Y[(z-y i ) a > 

i=i i=i 

where p and q do not have common roots. Then x and y are critical points of the Energy function 

i=n,j=n n i=n,j=m i=m,j=m ifi 

H(x,y)= ^2 VjViwjxi, Vjujxj)- ^2 ViO-jw(yi,Xj)+ a i ajw{y i , y 3 )-^ J Qjujyj) 

i<j=l i=l »J=1 i<j=l j=l 

(36) 

/ \ , (si ~ x if 
w(x u x 2 ) = In — 

P(x 1 )P(x 2 ) 
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In other words, the total charge at point Xi equals to the difference of multiplicities of the 
corresponding root in p and q. 

Proof: It is straightforward to check that 

H\ [p, q] (z) _ H\ \p, q] (z) 

— ies? 



^ p(z)q(z) q(z)p(z) 

where res z=Xi stands for the residue of a simple pole in the point The residue is zero since p 
and q satisfy (35). 

By direct calculation we get 

» - —Iff - + % S 1 - + 

which is a derivative dTL/dxi of the energy (36). Repeating similar calculation for we complete 
the proof. 

The following proposition gives examples of equilibrium configurations corresponding to several 
generic cases of the figure 1. 

Proposition 4 Let I — %\ < < ik+i be a strictly increasing sequence of nonegative integers 

and let Qi(z) be classical orthogonal polynomials satisfying the hypergeometric equation 

(L + \ i )Q i (z) = 0, L = P{z)^ 2 +U{z)j- z (37) 

where (up to a linear transformation of z) 

(i) P(z) = l U(z)=bz 

(ii) P(z) = -z 2 U{z) =bz TU 



Then polynomials p and q 

P (z) = P(z)^ k ^w[Q l M,QM,---Q ik (z),Q tk+1 ( z )l 

q(z) = P(z)X k - 1 »W[Q il (z),Q ia (z),...Q ik (z)] 

(i) deg(p) = n = ij - \k(k + 1), deg(g) = m = ^ fc =1 i, } - \k{k - 1) 

(ii) deg(p) = n = ij, deg(g) = m = £)J =1 ij 

satisfy the bilinear hypergeometric equation (35) with \ nm given in proposition 2. 

W[ipi( y z)...ip k (z)] = det \\d?pi(z) / dz^ \\ in (38) denotes the Wronskian determinant. 
To prove the proposition we need the following lemma by Crum [11] 



(38) 
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Lemma 4 Let L be a given second order Sturm- Liouville operator 

with a sufficiently smooth potential uq, and let {ipi, ...ip k } be its eigenfunctions corresponding to 
arbitrary fixed pairwise different eigenvalues {Ai, ...Xk}, i-e. ipi G ker(L + \), i = l..k. Then, for 
arbitrary ip G ker(L + A) the function 

satisfies the differential equation 

d 2 



dp 



+ u k ((j)) + X) ip = 



with 

d 2 

u k = u + 2— \nW[ipi...ipk\ 

Proof of proposition 4: 

Changing variables as in (18) and making a gauge transformation 

L — > L = vLv^ 1 , — In v 



d<j> 2VP 
we get a formally self-adjoint operator 

f d U \ ( d U \ d 2 

Lo "U + 27fJ \d^-^p)-w + UQ 

with eigenfunctions 

ipi = vQi, (L o + X i )ij} i = 0, i = 0,1,2,... 
According to the Crum lemma the function 

WbPh 0) , ^ 2 (0) ,-4>i k 0) ; VWi (0)] _ vp_ 

is an eigenfunction of 

d 2 d 2 d 2 

+ u k , u k = u + 2—logW{^ il ,i) i2 ,...i) ik }=u -2—\og(v k q) 



(39) 



with the eigenvalue Aj fc+1 . 
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Deriving (39) and the last equation we used the following properties of Wronskians 



/ dz 



in(n-l) 



W[vf U ...!//„] = .../„], W[h{z{<t>))-fn{z{<t>))\ = ) W[h{z)-fn{z)\ 

It then follows immediately that 



(L k + A ifc+ i) 



where 



= = 

(40) 



. dU _UdP 
dz P dz 



It is clear from the statement of the proposition that is independent of z. 

Finally, changing the independent variable back to z (18), <fi = <p(z) , we arrive to the bilinear 
hypergeometric operator in the rhs of (40). The degrees of p and q can then be evaluated from 
the highest powers of Wronskians (38). 

Example: Consider, for instance, the equilibrium configuration corresponding to the sequence 

7 = 2,4,6 (41) 

in the system with 

P = 1, U = -2z 

The eigenstates of the linear problem are Hermite polynomials H n (z) (see figure 1). Computing 
p and q, with the help of (38), we obtain 

p = 8l92z 3 (8z 6 - 12z 4 + 18^ 2 - 15), q = 32z(<iz 4 + 3 - Az 2 ) 

Polynomial p has a multiple root z = and this is a common root with polynomial q. 
Excluding common factors we have 

p = 256z 2 (8z 6 - l2z A + 18z 2 - 15), q = Az 4 - Az 2 + 3 

It can be verified without much difficulty that q and p do not have multiple roots, other than 
z = 0. Hence sequence (41) gives the following equilibrium distribution of charges ( interacting 
via logarithmic potentials ) , in the linear external field: One charge of the value v\ = +2 at z — 0, 
six charges of the value v 2 ..7 = +1 on the real line and four negative charges <7i.. 4 = —1 in the 
complex plane. 

The following proposition is an analog of proposition 4 for P(z) = z. 
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Proposition 5 Let I = i\ < %2---ik < h+i be a strictly increasing sequence of nonegative integers 
and k = mod 4. Let Qi(z) = L^f l \—bz) be Laguerre polynomials satisfying the hypergeometric 
equation 

(L + Xi)Qi(z) = 0, L = zj^ + bz-jt, b + 

Then polynomials p and q 



(42) 



p(z) = z^ k ^W[Q tl (z),Q l2 (z),...Q ik (z),Q tk+1 (z)l 
q(z) = zfr- 1 ^[Q il (z),Q i2 (z),...Q ih (z)] 

k+l ^ k 1 

deg(p) = n = ^^ij — ~^k(k + 1), deg(g) —m — ij — -^k(k — 1) 

3=1 3=1 

satisfy bilinear hypergeometric equation (35) with P(z) = z and U(z) = bz. 

Proof: repeats the proof of proposition 4, except that now k must be multiple of 4 in order to 
(42) be polynomials. 

VI. Rational and Trigonometric solutions of KP/KdV hier- 
archies, Evolution in Two Dimensions 

Another interesting set of examples are limits U — 0. They correspond to decreasing at infinity 
rational or periodic soliton solutions of the KP/KdV hierarchies (see equation (26) ). 
For instance, studying the case 

P = l, U = 

Bartman [4] provided an electrostatic interpretation for the Adler-Moser polynomials. 

Indeed, in this limit the bilinear hypergeometric equation becomes the recurrence relation for 
the Adler-Moser polynomials 

p"q - 2p'q + pq" = (43) 

which, as shown by Burchnall and Chaundy (who, by the author knowledge, first studied (43) in 
[7]), exhaust all polynomial solutions of (43). 

Note that, in difference from the generic cases shown in the figure 1, we have a set of polynomials 
depending continuously on k + 1 parameters: 

P — Qk+i, Q — 6k, @k = dk{z + ti, t 2 , ...tk) 

0fc = W[^i,...^fc], i>] = ^-i, il> = l,il>i = z, deg(6 k (z)) = (k + l)k/2 

with the second logarithmic derivatives of 8s being rational solutions of the KdV hierarchy. 

Thus, we have (at generic values of ti) equilibrium of k(k + l)/2 positive and \{k + l)(k + 2) 
negative free charges with positions in the complex plane continuously depending on U. 



17 



Let us turn now to the following problem: Find homogeneous polynomials p(X, Y, t), deg(p) = 
n , q(X, Y, t), deg(g) = m in two variables X, Y satisfying the equation 

f f q-f t P= (X 2 + Y 2 ) (qAp - 2(Vg, Vp) + P Aq) , (44) 

where 

A- — — — — 

and (, ) stands for the standard scalar product in C 2 . 
Factorizing p and q as 

n m 

p = Y\_(X sinifii — Y cos (pi), q = J^(Xsin6>j — YcosOi) 
i=i i=i 

we come to the following 

Proposition 6 The bilinear evolution (44) induces dynamics 



d0i 



(45) 



Proof: It is convenient to write (44) in the polar coordinates (X = r cos 0, Y = r sin 0) 

n m 

p = r n p = r n Y\ sin(0 - 0;), q — r m q = r m sin(0 - 6>j) 
i=i i=i 

_dp _dg ~<9 2 p <9g<9p <9 2 g _ 2 __ 

g ^ - ^ = q W " ^ + ^ + (n - m) M (46) 

Then, it can be verified that equation (46) corresponds to the case P(z) = —z 2 , U(z) = 
in the classification of figure 1. It must be remarked, however, that solutions p and q are not 
polynomial, but algebraic functions of z — exp(2i0): 

V = z ~ n/2 Yl je N x J 1/2 Pi x j = exp(2i0 j ), p = U j( , N ( z ~ x i) 

q = z -m/2 Yl j( _ M yj 1/2 q, Vj = exp(2i%), q = ]\ jeM (z - Vj ) 

Nevertheless, since p and q are of "almost" polynomial type, we get (45) by arguments similar to 
the proof of the proposition 2. More precisely, for this purpose it is rather more convenient to use 
(25), where we can replace q and p with "pure" polynomials p and q (and ip — | with |) : This 
substitution adds constants to coefficients in (25) and rhs of the equation acquires the common 
factor Z2( n ~ m ' ) exp(iJ2 jeN <f)j — iJ2 ieM dj). The lhs of (25) acquires the same factor, since the 
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quantity ("center of mass of the system") J^j^N^j ~~ X^eM^i does no ^ change with time. The 
later statement can be easily verified adding equations of motion for 0s and subtracting equations 
of motion for 9s in (46). Thus, the problem is reduced to the purely polynomial dynamics, which 
completes the proof. 

The flow (45) is a trajectory of two Sutherland systems in the absence of the external potentials. 
It is interesting that the equilibrium condition for (46) written in coordinates X, Y 

qAp - 2(Vg, Vp) + pAq = (47) 

has been studied in [6], [5] in connection with the Hadamard problem in the Minkowski space: 
Solutions of (47) (or fixed points of (45)) define differential operators possessing Huygens property 
in the Hadamard sense [13]. The angular parts p, q of solutions to (47) are periodic soliton solutions 
of the Korteveg-de Vries equation. 

The following proposition provides us with k + 1-parametric family of solutions to (47) describ- 
ing equilibrium configurations on the Coulomb charges (vortices) on the cylinder. 

Proposition 7 Let I — i\ < %2---ik < h+i be a strictly increasing sequence of nonnegative integers. 
Then 

P = r n W[il> h ,il> i2 , ...il> ik ,il> ik+1 ], q = r m W[il> h ,il> i2 ,...il> ik } 

where 

k k+l 

:= sin(i,-0 + U), W := det \ \dtpi/d4^\\, m = ^i,-,n = 

3=1 3=1 

satisfy (4?) 

Proof: repeats the proof of the proposition 4 for (46) =0, except that now we have a superpo- 
sition of the Thebyshev trigonometric polynomials sin(j0) and cos(j0) [20] instead of Qj. 

Thus, in contrast with the plane case we have k + l continuous parameters and k + l integers 
defining equilibrium configurations on the cylinder. Also, in difference from the plane distributions, 
the equilibrium is possible not only for consecutive triangle powers of the Adler-Moser polynomials, 
but for any two values of partitions n = deg(p) = Yl^lh an d m = deg(g) = Y^=ih- This is 
due to a different topology of the problem: roughly speaking, the charges on the cylinder have 
less "possibilities" to "escape" to infinity then on the plane. It must, however, be restated that 
numbers and values of charges depend on multiplicities and common factors of p and q. 

VII. Polylinear Evolution Equations and Related Hamilto- 
nian Systems 

As was mentioned in the introduction, (3) is, in fact, a special case of more general polylinear 
equation. The polylinear equation induces a polynomial dynamics which can be also embedded in 
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a Hamiltonian flow. However, this flow does not separate now in independent components. The 
Hamiltonians are of the Calogero-Moser type for several species of interacting particles. They are 
not generally related to the Coxeter reflection groups. 

Let us begin by introducing / species of particles with distinct charges Q := {Qi,i = 1../}, 
Qi 7^ Qj,i 3 ■ We define the /-linear differential operator operator 

itfbi, ...«K*) = p(z) (zLQi^mU^) +^ i<j QiQjq^(z)n l n^ j <in(z)) 

(48) 

+ \P'{Z) E L QU(Z) nU Qn(z) + U(Z) Ell QiQi(z) nU Qn(z) + A nil ft w 

acting on polynomials q ri e Vi = C H+1 ,i = I.. I 

Ni+rii i-l 

Qi(z)= II #i = I>j (49) 

j=Ni+l j=l 

For convenience, we now use unique numeration for roots of all polynomials. 

Similarly to the case I — 2, Qi — 1, Q 2 — —1 of the bilinear hypergeometric operator (f), 
the following proposition holds (We skip proofs below, since they repeat arguments of preceding 
sections) . 

Proposition 8 The polylinear operator n : V\ x V2 x ... x Vi — > Vj+i, Vj+i = C^ i=ini given 

by (48) with 

P(z) = A + Bz + Cz 2 , U(z)=a + bz, 

and 

\ m ,., ni = -(u'+ l -p" ) 

V ~ i=\ J 1=1 

induces dynamics 

^ = -2P(z t ) ^ - U(z t ) - 9lp\ Zi ) (50) 

by action 

i=l n^i 

on 

Under conditions mentioned before, (50) describes motion/equilibrium of n-y charges Q±, n 2 charges 
Q2 etc. 

To avoid confusion, we note that in (50) and up to the end of this section, the summation 
indexes go from 1 to the total number of roots Ej=i n j an d to each root Zj we assign Qj which is 
equal to the charge of the corresponding polynomial. 

It is natural to ask the following important question: May dynamics (50) be embedded in a 
Hamiltonian flow? We address it in the following 
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Theorem 3 (50) is a trajectory of the Hamiltonian system 



H 




where 



W{ Zl ,z 2 ) 



P( Zl ) + P(z 2 ) 
(zi - z 2 f 



U Qj (z) = U(z) + QjP f (z)/2 



System (51) is of the Calogero(Sutherland)-Moser type of I species of particles with masses Qi. 
In this picture, the two-body potentials are of the similar form, but having different amplitudes 
QjQk{Qj + Qk) (for interaction within each of the species of particles and between the species 
respectively). They are translation invariant in Newtonian coordinates (18) if (up to a linear 
transformation) P(z) = 1 or z 2 . We remind that in the case I — 2, Q\ — 1, Q 2 — —1 the interaction 
between two different species vanishes, leading to separation of the Hamiltonians, while for I — 1 
we obtain identical Calogero-Moser particles. Both above cases are related to A/BC/D root 
systems. 

It is seen without much difficulty that for generic /, Q, system (51) is not related to any 
Coxeter reflection group. Although , by the author knowledge, quantum models related to different 
deformations of the Coxeter root systems were considered in earlier works ( eg [6], [9] ), (51) has 
not appeared in the literature. We do not attempt to address the question of integrability of (51) 
in this paper, leaving it for future studies. 

VIII General Bilinear Dynamics, Evidences of Integrability 

Let us concentrate on the two-component case I — 2, Q — 1, — A of general bilinear hypergeo metric 
operator (23), restricting ourselves with the rational case 



of n and m particles of two different types. According to proposition 8 and theorem 3, (52) is a 
corollary of the evolution equation 



P(z) = i,U (z) = iujz, i := V^l 



Such a choice of the coefficients leads to dynamical system 




(52) 




iq-f- - iAp-f = p"q - 2Ap'q' + A 2 pq" + tuz (p'q - Apq') + u(Am - n)pq 



(53) 
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for polynomials p and q (22), and is a trajectory of a system with the Hamiltonian 

jr -i£((S , H^£((*) ,+ ^) +vfe ' ) (54) 



V 2 A(A- 1) 2A 3 

x ' y ~ fa - Xk)2 + - y^ 2 khtM to - y ^ 2 

Remark: Although, similarly to the case A = 1, (53) can be written in the form of a time 
dependent Schrodinger equation (25), (changing t to it) with the "■?/>" and "r" functions given by 

^ = j,/ g A )T = f A(M 

the dynamics of poles of the potential U = 2(logg)" cannot be embedded in a Hamiltonian flow 
uncoupled from the dynamics of zeros of p. This is why system (52) cannot be connected with 
solutions of the KP hierarchy. 

Although the Hamiltonian system (54) unlikely be integrable for arbitrary initial conditions 
and A ^ 0, ±1, we find that its trajectories (52) (defined by the polynomial evolution (53)) are 
integrable. 

Conjecture 1 System (52) is completely integrable for arbitrary real A and u in the sense that 
there exist 2(n + m) — 1 functionally independent integrals of motions Ij, which are real rational 
functions ofx,y, i.e. Ij = Ij(x u x\, ...x n , x* n , y ± , y\, .., y n , y*) = Ij(x\, x u .., x* n , x n , y{, y u .., y*, y n ), 
j = \..2{m + n) - 1 

We devote the rest of this section to examples in favor of this conjecture. 
We take the case A = 1 as the first example: According to corollary 1 the equations of motion 
can be reduced to the Lax form (see eg [19]) 

d d 
i ~dt Lx = ][Lx ' Ax ^ + U}Lx ' i ~dt Ly = ]yLy ' Ayi + wLy 

where L y , A y and L y , A y are matrices of dimensions n x n,m x m respectively 

(L.hk = i(i% + vx,) S kj + i=4t j,keN 

(Ly)jk = \ {it + <"Vi) ^ + J,keM 



(55) 



Substituting (24) to (55) we eliminate velocities dx/dt, dy/dt, getting the Lax matrices L x and L y 
for (24), depending on the coordinates only. It is easy to see that the absolute values of squares 
of traces 

h = (TrL)(TrL)*, I 2 = (TrL 2 )(TrL 2 )*, .., 7 2(n+ro) _i = (TrL 2 ^ +m ^~ 1 )(TrL 2 ^ +m ^ 1 )* (56) 
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Figure 2: Examples of trajectories of the two component system consisting of n = 6 unit charges 
and m = 1 charge —A, where A = 1.213579. Every curve shows an individual trajectory of each 
charge in its own coordinates ($lxj, ^sxj),j = 1..6 or (dtyj, $syj),j = 1. The charge of value —A is 
depicted by the gray solid line. The motion shown on the left figure has period An/oo = 2T, while 
the period on the right figure is equal to T = 2tv/uj. The motion on the left is depicted within the 
time interval t = [0, T], which is a half-period for such initial conditions. In this case trajectories 
of several charges coincide interchanging each half-period T. 
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of the (m + n) x (m + n) matrix 



L x 
Z„ 



are real rational integrals of motion. They are homogeneous functions in x, y, u (with x, y and ui 
having weights —1, —1,2 respectively). The functional independence of (56) can be easily proven 
by considering them as polynomials in uj with functionally independent highest symbols 



h = w * > x + > ^ y > >*)* + > (y*r ]+., k= 1..2(n + m) - 1 




Let us now turn to the general system A ^ 0, ±1. Although, in this case, our arguments in favor of 
inetgrability of (52) stem mainly from numerical studies, we would like to mention some analytic 
results: 

The Hamiltonian (54) admits total separation of variables in low dimensions n + m < 4. 
Namely, separating motion of the center of mass, we obtain one or two dimensional problem, 
admitting (in the latter case) further separation of variables in the polar coordinates. 

Another (less trivial) example is the system with even number of unit charges n = 21 and a 
single particle of the second type m = 1 having an arbitrary charge —A. The system is subject to 
symmetric initial conditions: 

Xj (p)=x j+l (0), j/i(0) = 0, j = l..l 

It is seen without much difficulty that due to this Z 2 symmetry, the above conditions hold for 
any t. Taking this fact into account, we may reduce (52) by this symmetry keeping only variables 
Xj,j = Changing the variables as Xj = J~z] we arrive to the following equations of motion 



i d 4l = -Y^-+u (57) 



k+3 

which corresponds to the BC n rational case of figure 1. The integrability of (57) is then provided 
by arguments used for the study of linear dynamics. The rational integrals of motion for (57) can 
be found using the Lax representation for the Calogero-Moser system of the BC„ type. 

One can also prove the periodicity of small nonsymmetric deviations from the symmetric 
trajectories as linear perturbations around (57). We do not perform this analysis here, since it 
requires cumbersome calculations. 

Finally, numerical simulations show that for any initial conditions and real A trajectories of 
(52) turn out to be periodic, which shows existence of 2(n+m) — 1 independent rational integrals of 
motion. The period of motion is an integer multiple the period of the "free" oscillator T = 2-k/lu, 
with an integer factor depending on the initial conditions. 

Typical examples of trajectories for generic initial conditions and A are shown on figure 2. 



24 



IX. Conclusions and Open Questions 



The results at which we have now arrived may be summed up as follows: The bilinear hyperge- 
ometric operator (1) induces dynamics (2), which may be embedded in a Hamiltonian flow. In 
the case A = 1 this flow is generated by a sum of two independent Calogero(Sutherland)-Moser 
Hamiltonians (theorem 2) with (generally) different forms of external potentials. This allows us 
to integrate (24) by the Lax method. The fixed points of the bilinear evolution correspond to 
equilibrium distributions of different species of point vortices on the plane or cylinder. They may 
be obtained (again, for A = 1) by a finite number of Darboux transformations from the eigenstates 
of associated linear problems (propositions 4, 5, 7). The dynamical system (52) of two species of 
interacting points in an external field is conjectured to be completely integrable for arbitrary real 
A and uj. 

Let us now mention some open questions 

The main problem, of course, is to prove integrability of (52) for arbitrary A (conjecture 1). It 
might be done by using two approaches: The first approach is to find a Lax representation for (52). 
The Lax matrices for arbitrary A could have a complicated structure, being rational functions of 
x, y and u of greater homogeneity degree in comparison with the Calogero-Moser case. It makes 
hard to find them using straightforward computational approaches. 

Another method is to try to linearize (53). Although, as was mentioned before, such a lineari- 
sation, connected with the KP equation, was possible for P = 1 and A = 1, we cannot apply a 
similar scheme in general case. 

Another set of questions is connected with solutions of bilinear hypergeometric equation: 

In 1929 [7], Burchenall and Chaundy studied the following question: What conditions must be 
satisfied by two polynomials p(z), q(z) in order that the indefinite integrals 



may be rational, provided p and q do not have multiple and common roots? 

They found that the integrals are rational if p — 9 i: q — 9 i+ i are (now know as the Adler-Moser) 
polynomials satisfying (43). 

From the other hand, we know that any two polynomial solutions of the ordinary hypergeo- 
metric equations p = Q n , q = Q m are orthogonal with the measure u(z) 



Since the two above integrals are related to particular forms of (4), it is natural to ask the following 
question: what integration condition may be imposed on two polynomials p and q in order for 
they satisfy a nondegenarate bilinear hypergeometric equation (4)? 

In the same work [7], Burcnall and Chaundy have shown that any polynomial solutions of 
(43) may be obtained by a finite number of Darboux transformation from the kernel of the "free" 
differential operator d 2 /dz 2 . 
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In this paper we have proved proposition 4, stating that polynomials obtained from the eigen- 
functions of the ordinary hypergeometric equation by a finite number of Darboux transformations 
are solutions of (4). By analogy with [7] it is natural to state the following 

Conjecture 2 Any polynomial solutions to bilinear hypergeometric equations of propositions 4, 5 
are (38) and (42) respectively. 

Concluding the article we would like to mention briefly possible multi-dimensional general- 
izations of (4). A particular generalization was constructed in section VI for the homogeneous 
polynomials in two variables (47). A similar construction [5] related to the classical special func- 
tions in many dimensions [14], is connected with the quantum Calogero-Moser systems on the 
Coxeter root systems (and their deformations [6] , [9]). In this context, it would be interesting to 
find a proper analog of (4) in many dimensions. 
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